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Abstract 



We provide a theoretical foundation for the statistical approach for computing the absorp- 
tion properties of particles in the Rayleigh domain. We present a general method based 
on the Discrete Dipole Approximation (DDA) to compute the absorption and scattering 
properties of particles in the Rayleigh domain. The method allows to separate the geomet- 
rical aspects of a particle from its material properties. Doing the computation of the optical 
properties of a particle once, provides them for any set of refractive indices, wavelengths 
and orientations. This allows for fast computations of e.g. absorption spectra of arbitrar- 
ily shaped particles. Other practical applications of the method are in the interpretation of 
atmospheric and radar measurements as well as computations of the scattering matrix of 
small particles as a function of the scattering angle. In the statistical approach, the optical 
properties of irregularly shaped particles are represented by the average properties of an 
ensemble of particles with simple shapes. We show that the absorption cross section of an 
ensemble of arbitrarily shaped particles with arbitrary orientations can always be uniquely 
represented by the average absorption cross section of an ensemble of spheroidal particles 
with the same composition and fixed orientation. This proves for the first time that the 
statistical approach is generally viable in the Rayleigh domain. 
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1 Introduction 



The absorption and scattering properties of small particles are very important in 
both astronomical and atmospheric remote sensing applications. The interaction of 
light with particles much smaller than the wavelength of radiation has been stud- 
ied first by Lord Rayleigh [1] who explained from basic physical principles the 
color and polarization of the light from the sky. When studying the detailed spec- 
tral properties of the interaction of particles with light, we have to consider the 
effects of particle size, shape and composition. In this paper we discuss the effects 
of shape and composition of homogeneous particles, while taking the particle sizes 
to be in the Rayleigh domain, i.e. much smaller than the wavelength both inside 
and outside the particle. 

When considering the possibilities for computing the optical (i.e. the absorption and 
scattering) properties of small particles there are two extreme approaches one can 
take. The first approach is to assume the particles are homogeneous in composition 
and spherical in shape, which allows us to perform fast and simple computations 
of its interaction with light using Mie theory [2]. The second extreme is to make a 
model of the particle in an exact way, and perform numerical calculations to obtain 
its optical properties. This can be done using, for example, the so-called DDA (Dis- 
crete Dipole Approximation; see e.g. [3]) or the T-matrix method (see e.g. [4,5]). 
The first approach is fast, and provides insight into the physics and effects that play 
a role in the interaction of light with small particles (see e.g. [6]). However, due 
to the perfect symmetry of homogeneous spheres, resonance effects may occur, for 
example, at particular values of the refractive index, that are not seen in realisti- 
cally shaped natural particles. This limits the applicability of this approach. The 
second approach allows us to reproduce details in the observed properties of irreg- 
ular particles. The main drawback is, however, that the computational demand of 
most numerical techniques available to compute the optical properties of realisti- 
cally shaped particles, is high. If we wish to consider a large collection of various 
particle compositions, sizes or wavelengths, one can resort to a third method, the 
statistical approach. In the statistical approach one simulates the average optical 
properties of an ensemble of irregularly shaped particles by the average proper- 
ties of an ensemble of particles with simple shapes. These simple shapes guarantee 
computations of the optical properties to be relatively fast. In addition, by choosing 
a broad distribution of simple shapes, we can get rid of the resonance effects which 
can occur when using homogeneous spheres. The statistical approach is proven 
useful for, for example, computing absorption spectra of small forsterite grains [7], 
and for calculating the degree of linear polarization for small quartz particles [8]. 
Some of the successes and limitations of the statistical approach are discussed in 
refs. [7,8,9,10,11,12]. 

In this paper we employ an analytical method based on the DDA to compute the op- 
tical properties of particles in the Rayleigh domain in an efficient way. This method 
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is then used to show that the average absorption cross section of an ensemble of 
particles with arbitrary shapes and orientations is identical to that of an ensemble 
of particles with spheroidal shapes in a fixed orientation and with the same compo- 
sition. We thus provide an analytical basis for the use of the statistical approach in 
computing the absorption cross sections of particles in the Rayleigh domain. 

In Sect. 2 we outline the method which is summarized for practical purposes in 
Sect. 2.4. In Sect. 3 a few examples of applications to complex particles shapes 
are provided. A discussion of the implications for the applicability of the statistical 
approach is given in Sect. 4. In this section we will also provide suggestions on how 
to use the derived shape distribution of spheroids to calculate the optical properties 
of particles with sizes outside the Rayleigh domain. 



2 The absorption properties of very small particles 



In this section we will outline the method to obtain the optical properties of arbitrar- 
ily shaped and arbitrarily oriented particles small compared to the wavelength. We 
do this using a solution of the DDA. The general equations for the optical properties 
of particles in the Rayleigh domain are given in Sect. 2.1. The method to compute 
the optical properties for arbitrarily shaped particles is outlined in Sect. 2.2 and 
summarized in Sect. 2.4. In Sect. 2.3 we introduce the distribution of form-factors. 
Using this distribution we prove the following fundamental theorem: 

Theorem 1 The absorption cross section of an arbitrarily shaped and arbitrarily 
oriented, homogeneous particle in the Rayleigh domain, or an ensemble of such 
particles with various shapes and orientations, equals the average cross section of 
an ensemble of spheroidal particles in a fixed orientation with the same composi- 
tion and a shape distribution that is independent of the composition of the particle. 

This theorem proves the validity of the statistical approach for calculations of the 
absorption cross sections for particles in the Rayleigh domain. It is also an encour- 
aging result for applications of the statistical approach for other purposes. 



2.1 General equations in the Rayleigh domain 



In the Rayleigh domain a particle in a given orientation interacts with incident light 
as a single dipole with dipole moment p given by 

p = orE inc . (1) 
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In this equation a is the 3x3 polarizability tensor, and Ei nc is the incoming electric 
field. From the polarizability tensor it is possible to obtain all scattering and absorp- 
tion properties of the particle (see e.g. [13]). The absorption cross section of such 
a particle when the incident field is applied along the x, y or z axis of a Cartesian 
coordinate system (indicated by the symbol fx = 1,2, 3) is given in ref. [14] 

C^=kJm(a^), (2) 



and the scattering cross section by 




where a^ v) is the ji, v-th component of the polarizability tensor (/i, v = 1,2, 3), and 
k = 2n/A, with A the wavelength of incident radiation. Note that when a particle 
is in the Rayleigh domain the electric field can be taken constant over the particle 
volume. This implies that the direction of wave propagation of the incident light is 
not important for the interaction of the light with the particle, but only the direc- 
tion of the electric field. The cross sections averaged over all particle orientations 
(denoted by (...)) are simply given by the average over the three axes which results 
in 



(C abs ) = 3 Im 
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(4) 



and 



\a 



G"v)| 



(5) 



For only a few particle shapes simple analytical equations are available to compute 
the polarizability tensor. These are the shapes of homogeneous, layered and hol- 
low spheres, and homogeneous, layered and hollow ellipsoids. For a homogeneous 
sphere a is given in ref. [13] 



m 2 - 1 
a = 3V^— I 3 



(6) 



with m the complex refractive index of the particle material, V the material vol- 
ume of the particle, and I 3 the 3 x 3 identity matrix. The polarizability tensor of a 
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homogeneous ellipsoidal particle is given by ref. [13] 



a = V 



m L -\ 
1+Li(m 2 -1) 








m 2 -l 
1+L 2 (m 2 -1) 








m 2 -l 
l+L 3 (m 2 -\) J 



(7) 



where the L, are geometrical form-factors with values between and 1 depending 
on the shape of the ellipsoid. 

For an ellipsoid of revolution (a spheroid) with axes a and b we have 



U = 



1 1 , < 1+e 
-1 + — In 
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and L 2 = L 3 = (1 -L0/2. Here b is the rotation axis, and a/bis the aspect ratio. For 
prolate spheroids a/b < 1 and < L x < 1/3, while for oblate spheroids a/b > 1 
and 1/3 < Li < 1. For homogeneous spheres L x = L 2 = L 3 = 1/3. 

The equations for the polarizability of hollow spheres or ellipsoids can be found 
in ref. [6]. For a hollow sphere having a volume fraction / occupied by a central 
vacuum inclusion, the polarizability tensor is given by 

(m 2 - l)(2m 2 + 1) 

a ~ V (m 2 + 2)(2m 2 + 1) - 2(m 2 - Iff 3 ' ( ) 



where V = |^r 3 (l - /) is the material volume of a hollow sphere with outer radius 
r. 

For convenience we introduce the dimensionless polarizability per unit material 
volume of a particle a = a/V. 



2.2 The Discrete Dipole Approximation ( DDA ) 



In the Discrete Dipole Approximation (DDA) a particle is represented by a col- 
lection of interacting dipoles. This idea was proposed by Purcell & Pennypacker 
[15], who derived the general equations for the interaction of the dipoles with the 
incident light and with each other. The approach is equivalent to a discretization of 
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the particle volume and the assumption that each volume element interacts with the 
radiation field as a single dipole j with polarizability fij. For a theoretical derivation 
of this approach from the Maxwell equations see e.g. [16,17]. For inhomogeneous 
particles the polarizability can vary throughout the particle depending on the ma- 
terial properties of the local volume element. Also the polarizability tensor can be 
anisotropic. In this paper we consider only homogeneous particles, i.e. fij = fi for 
all j. In addition we consider only the case when the polarizability is isotropic in 
which case the polarizability tensor can be written as /3I 3 , with /? the scalar polar- 
izability. For the DDA to be valid, the volume elements represented by the dipoles 
have to be in the Rayleigh domain. 

To compute the interaction of N dipoles with the incident field and with each other 
we have to solve a set of /V inhomogeneous linear vector equations, namely 

N 

J^iAij +/3'%)Pj = E incJ (i = 1, AO. (10) 



Here the polarization vector at the position of dipole j, P ; , is a column vector with 
3 elements. E inCji is a column vector representing the incoming electric field at the 
position of dipole i. The Ay are 3 x 3 matrices determined by the positions of the 
dipoles in the particle and thus by the particle geometry. The general equations for 
the Ay are given in ref. [3]. We assume in this paper that the entire particle is in the 
Rayleigh domain. This implies that the electric field is constant, so we are in the 
limit of k — » 0. Then 



h - 3f,;,f, 



'J"-'J 



Anrl 



0, 



i = J- 



(11) 



Here r ; y is the unit vector pointing from dipole i to dipole j, and r (; is the distance 
between these two dipoles. The 3x3 matrix r, 7 r iy is defined as the product of r !; 
as a column vector and r i; as a row vector, i.e. 



~2 

rr. 

IJ,X 



r ij,x r ij,y 



r ij,x r ij,z 



r ij,x r ij,y 



r ij,x r ij,z 



r >j>y r U,z 



r 2 . 



(12) 



where r ij>x , r^ y and r^ z are the x, y and z components of the unit vector f iy -, respec- 
tively. 



6 



Following [3] we define the 3N dimensional column vectors Ei nc = (Ei nC; i, E; nCj2 , 
Einc,iv) andP = (Px,P 2 , ...,P;v) and the 3N x 3./VmatrixA such that A 3(i _ 1 ) +jUi3 ( / _ 1)+v = 
(Ay)^ v , {/, j = I...N; n, v = 1, 2, 3}. In this way Eq. (10) can be written as a single 
matrix equation 

(A+fi-% N )P = E inc . (13) 

Note that in this equation the matrix A only depends on the particle geometry, while 
the matrix /T 1 ^ only depends on the dielectric properties of the particle. 

Eq. (13) can be solved for P using, for example, a conjugate gradient method. In 
the literature considerable effort has been reported to increase the efficiency of so- 
lution methods for Eq. (13) (see e.g. [18,19]). In contrast to the general equations 
for A, in the Rayleigh domain this matrix is real and symmetric. This implies that 
the eigenvalues of A are real, and that the eigenvectors form an orthonormal basis. 
As we will detail below, once the eigenvalues and eigenvectors of A have been de- 
termined, the solution of Eq. (13) can be obtained for arbitrary values of /3, and thus 
for any particle material. This method is similar to that presented by Markel et al. 
[20]. 

Using the eigenvalues and eigenvectors, the matrix A can be diagonalized by writ- 
ing 

A = UDlf, (14) 

where a superscript T is used to denote the transpose of a matrix. The columns 
of U are the 3N real valued eigenvectors of A, and the elements of the diagonal 
matrix D are the corresponding eigenvalues of A. Since the eigenvectors of A are 
orthonormal the transpose of U equals its inverse, i.e. U T = U _1 . The eigenvalues 
and eigenvectors of a real symmetric matrix can be found relatively easy using 
numerical techniques. 

The matrix on the left hand side of Eq. (13) can be rewritten after diagonalization 
since 

(A + /T 1 I 3iV ) = (UDU r +P~% N ) = U(D +/rV)U r (15) 

Using this, the solution of Eq. (13) can be written as 

P = U(D+^- 1 I 3i vr 1 U r E mc . (16) 

Since (D + fi~ l I 3N ) is a diagonal matrix, its inverse is readily computed. Once U 
and D have been calculated, the solution can be easily computed for every value 
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of /3. Since U and D only depend on the geometry of the particle, we only have to 
calculate them once in order to find the solution for arbitrary values of J3, and thus 
for any particle material. 

In the Rayleigh domain, the incident electric field constant throughout the particle. 
This means that, at each dipole j, E incj = E inc . The total dipole moment of the 
particle then is 

N 

p = ^P y = <*E mc (17) 



Combining this equation with Eq. (16) gives us the polarizability tensor a with 
elements 



3N 



w 



G«v) 



(18) 



where the Aj are the 3N eigenvalues of A, and the wj v) are given by 

N 

W T = Z U 7 -,3( ( --l) + ,U 7 -,3(,-l) + v • (19) 



Lk=l 



Note that since the matrix U is real valued, the wj are real valued and positive 
since 



N 



W 



V (=i 



(20) 



Using Eqs. (2) and (18), the absorption cross section when the field is applied along 
the n axis becomes 



( 3N 



(21) 



Also all the scattering properties, i.e. the scattering cross section, matrix, and depo- 
larization factors, of the particle can be obtained from the polarizability tensor. The 
scattering cross section when the field is applied along the jx axis can be computed 
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using Eqs. (3) and (18) and becomes 



74 3 
^sca — ^_ / , 



6n 



v=l 



3N 

I 



w 



,0"v) 



(22) 



The orientation averaged cross sections can also be obtained from Eq. 18 using 
Eqs. (4) and (5). In the following we focus on the absorption cross section. 



2.3 The distribution of form-factors 



In the previous sections we have shown how to separate the geometrical and di- 
electrical properties when computing the absorption properties of arbitrarily shaped 
particles in the Rayleigh domain. In this section we show how the particle geome- 
try can be uniquely described by a distribution of form-factors. Also, we will show 
how this distribution can be used to construct a shape distribution of spheroidal 
particles in a fixed orientation and with the same composition which has identical 
absorption properties. 

We take the polarizability of a single volume element to be the Clausius-Mossotti 
polarizability of a cubic volume element with size d [15], i.e. 

= 3d 3 . (23) 

m l + 2 



When the size of the volume elements is sufficiently small, this is the polarizability 
that follows directly from the Maxwell equations (see e.g. [16,17]). 



The equation for the absorption cross section (Eq. 21) then reduces to 

m 2 -l 



-abs 



3/V (uu) 



7=1 



kV Im 



1 + Lj (m 2 - 1) 



(24) 



with V = Nd 3 the total material volume of the particle and the form-factor Lj 
defined by 

Lj = i + Ajd \ (25) 



Comparing Eq. (24) with Eq. (7) we see that Eq. (24) is equal to the equation for the 
average absorption cross section of an ensemble of oriented ellipsoids with various 
shapes. We have thus proven theorem 1. 
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The result given above implies that we can represent the absorption properties of an 
arbitrarily shaped and arbitrarily oriented particle uniquely by a distribution of the 
form-factors, P(L). But in order to represent an ellipsoid with a given aspect ratio, 
the Lj have to be between and 1. This implies that the eigenvalues, Aj, have to be 
between -l/(3d 3 ) and 2/(3d 3 ) for arbitrary shapes. For all particle shapes that we 
considered this was the case. From the w { j MM) (given by Eq. 20) and the Lj (given by 
Eq. 25) it is straightforward to construct the distribution of form-factors. 

In practice we often meet ensembles of particles with different shapes and orienta- 
tions. In the Rayleigh domain, averaging over particle orientations can be done by 
averaging over three perpendicular orientations. In view of Eq. (24) it is clear that 
the shape and orientation averaged polarizability (and thus the average absorption 
cross section) of such ensembles can be obtained from the ensemble averaged form- 
factor distribution. This allows us to construct average form-factor distributions of 
different classes of particle shapes, resulting in a unique distribution of spheroidal 
particles to represent the absorption properties of each class. 

2.4 Summary of the procedure 

The method described above provides a practical way of computing the optical 
properties of arbitrarily shaped and arbitrarily oriented particles in the Rayleigh 
domain. It separates the geometrical and material properties of the particle. This 
implies that for a given particle geometry the computation has to be done only 
once to obtain the absorption and scattering properties for arbitrary values of the 
refractive index and wavelength. 

The procedure to obtain the optical properties of an arbitrarily shaped and arbitrar- 
ily oriented particle in the Rayleigh domain is outlined as follows. 

(1) Construct the matrix A using Eq. (11) and A3(,-_i) +/i> 3( / -_i) +v = (Ay)^. 

(2) Compute the eigenvalues, Aj, and the matrix of eigenvectors, U, of the real 
symmetric matrix A. 

(3) Compute the wf v) using Eq. (19). 

(4) Compute the polarizability tensor a from Eq. (18). 

(5) From the polarizability tensor the absorption and scattering cross sections for 
the particle in a single orientation are obtained from Eqs. (2) and (3). Also, 
the scattering matrix can be obtained from the polarizability tensor (see [13]) 

(6) The absorption and scattering cross sections for an ensemble of particles with 
random orientations can be obtained from Eqs. (4) and (5). The scattering 
matrix for an ensemble of randomly oriented particles is given in ref. [13]. 

The computationally demanding steps (1), (2) and (3) have to be done only once 
for a given particle geometry and structure. Steps (4), (5) and (6) then provide the 
optical properties for arbitrary orientation, wavelength, and particle material. 



10 




Sphere 



Spheroid (a/b=0.5) 



H 1 HH 1 1 1 H — I 1 1 1 1 1 1 1 H 




J 1 1 1 1 L 



Spheroid (a/b =2.0) 



111 



\ ^-^^ 1 , , , 



0.2 



0.4 0.6 
L 



0.8 



Fig. 1. The form-factor distributions as calculated for an array of 3000 dipoles representing 
a sphere, a prolate spheroid with a/b - 0.5 and an oblate spheroid with a/b = 2.0 averaged 
over all particle orientations. The expected positions of the peaks in the distribution are 
indicated by the dashed lines. 

The distribution of form-factors is computed from the Aj and the wj^ by using 
Eq. (25). The corresponding distribution of spheroids is obtained from the form- 
factor distribution using Eq. (8). In Sect. 4.4 we outline how the distribution of 
spheroidal particles that is obtained might be used to compute the optical properties 
of particles outside the Rayleigh domain. 



3 Application to various shapes 



3. 1 Homogeneous spheres and spheroids 



Since the average polarizability of an ensemble of homogeneous spheroids in ran- 
dom orientation consists of the sum of two polarizabilities with different values of 
L (see Eqs. 7 and 8), the ensemble averaged form-factor distribution displays two 
distinct peaks, one at the value of L determined by Eq. (8) and one at (1 - L)/2 
(with double intensity). 

By way of example we computed the form-factor distributions of a homogeneous 
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Fig. 2. Contour plots of the imaginary part of the polarizability per unit volume as func- 
tions of the real and imaginary part of the refractive index. The left panel is for a single 
homogeneous sphere and the right panel for an ensemble of randomly oriented prolate ho- 
mogeneous spheroids with aspect ratios a/b - 0.5. 

sphere, a prolate spheroid with aspect ratio a/b = 0.5 and an oblate spheroid with 
a/b = 2.0, averaged over all particle orientations using the DDA method (see 
Fig. 1). The material volume of each particle is discretized by using 3000 dipoles. 
We also indicate in this figure the positions of the peaks as determined by Eq. (8). 
We see that, in general, the positions of the peaks in the calculated distributions 
match the expected values fairly well. The differences are mainly caused by the 
fact that by using a discrete array of cubic dipoles it is not possible to represent the 
shape of a spheroid exactly. 

The refractive index of natural materials is a function of wavelength. The refractive 
index of most materials exhibits resonances in the infrared part of the spectrum. 
These resonances lead to corresponding resonances in the absorption cross section 
which can be detected, for example, in thermal emission spectra. The wavelength 
positions of the features are sensitive to the particle shape and can be used to deter- 
mine the shape and composition of the emitting grains. Therefore, it is important 
to study the absorption cross section as a function of the refractive index and wave- 
length. 

In the left panel of Fig. 2 we show C a b s /(W) or, in other words, the imaginary part 
of the polarizability per unit volume, of a single homogeneous sphere as a function 
of the real and imaginary part of the refractive index as computed from Eq. (6). Note 
the sharp increase of the imaginary part of the polarizability when m —> i V2~ typical 
for homogeneous spherical particles (cf. Eq. 6). When considering the absorption 
cross section as a function of wavelength near a resonance, this leads to a sharp peak 
in the absorption spectrum. The imaginary part of the polarizability per unit volume 
of an ensemble of randomly oriented prolate spheroids with aspect ratio a/b = 0.5 
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is shown in the right panel of Fig. 2. In this contour plot we see two sharp maxima, 
caused by the two values of L determining the polarizability. In general, we expect 
particles with a few isolated form-factors to show strong resonances for each form- 
factor. As will be shown below, a full distribution leads to a smoother absorption 
spectrum. 



3.2 Gaussian random sphere 

In order to model the optical properties of irregularly shaped particles we have to 
employ a model for the shape of the particles. A successful model that is frequently 
used in light scattering theory is that of a Gaussian random sphere [21]. The ba- 
sis of this shape model is a homogeneous sphere of which the surface is distorted 
according to a Gaussian random distribution. The distortion of the surface is pa- 
rameterized by two shape variables, the standard deviation, cr, of the distance to the 
center, and the average correlation angle, Y. The value of Y determines the num- 
ber of hills and valleys on the surface within a solid angle, while cr determines the 
height of these hills and valleys. For details see ref. [21]. Using different shape pa- 
rameters we can create particles with varying degrees of irregularity. Volten et al. 
[22] successfully used Gaussian random spheres to compute the scattering matrices 
of mineral particles as functions of the scattering angle. They found that in order 
to reproduce the experimentally determined scattering behavior they had to employ 
very large values of cr. 

In the top panel of Fig. 3 we show the form-factor distribution of various ensem- 
bles of randomly oriented Gaussian random spheres with different values of cr. In 
all cases we fix the value of Y to 10°. For each value of cr we construct an ensemble 
of 10 different Gaussian random spheres, each with the same value for the shape 
parameter but different seeds of the random number generator. Above the two pan- 
els of Fig. 3 we show some typical examples from the ensembles. The volume of 
each particle is discretized using 3000 dipoles in order to construct the orienta- 
tion and ensemble averaged form-factor distribution. The form-factor distribution 
of each ensemble is used to compute its average polarizability per unit volume as 
a function of the refractive index from Eq. (24). The imaginary part of the average 
polarizability per unit volume is shown by contours in the bottom panel of Fig. 3. 

Fig. 3 shows that for small deformations of the surface of the sphere (most left 
panel), the form-factor distribution is dominated by one peak at L « 1/3 and the 
absorption properties are rather similar to those of spherical particles (see the left 
panel of Fig. 2). For larger values of cr the distribution of form-factors gets broader. 
If cr = 0.5 or 0.7 we see that for a fixed value of Re(m) the polarizability is a 
smooth, continuously increasing function of Im(m), as expected in the absence of 
resonance effects caused by the geometry of the particle. 
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Fig. 3. The orientation and ensemble averaged form-factor distribution (top panels) and 
the imaginary part of the average polarizability per unit volume as a function of the real 
and imaginary part of the refractive index (bottom panels) of the ensembles of randomly 
oriented Gaussian random spheres. Above the figures we show typical examples of particles 
from the various ensembles. From left to right the relative standard deviation of the distance 
to the center is 0.1, 0.3, 0.5 and 0.7, respectively. The correlation angle T = 10°, and is 
the same for all four shapes. For each shape we considered an ensemble of 10 different 
particles with the same parameters but different seeds of the random number generator. We 
average over these different particles in order to reduce the effects of random variations of 
the particle shapes on the final results. 



In ref. [23] a method is proposed to compute the scattering and absorption prop- 
erties of Gaussian random spheres in the Rayleigh domain. However, the method 
used in that paper is an approximation using only a single ellipsoid as a best fit 
to the particle shape. From the discussion given above it is clear that this can pro- 
vide only a very rough approximation to the real absorption properties of Gaussian 
random spheres. 
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3.3 Aggregates with various fractal dimensions 



In environments where dust particles grow by aggregation of smaller constituents, 
the resulting grains might be modeled as fractal aggregates. The fractal dimensions 
of these aggregates are a measure for the compactness or fluffiness of the particle 
and depend on the conditions in the environment. A large part of the interplanetary 
dust particles collected in the Earth's atmosphere are aggregates of small particles. 
Therefore, the optical properties of particle aggregates are very important in astro- 
physics and planetary physics. Fogel & Leung [24] computed emission and extinc- 
tion spectra of fractal aggregates with various fractal dimensions. They found that 
the mass absorption and extinction coefficients of fractal aggregates are on aver- 
age higher than those of volume equivalent spheres. Mackowski [14] developed a 
method to obtain the optical properties of sphere clusters in the Rayleigh domain. 
For computations of various values of the refractive index, which are needed, for 
example, for computing absorption spectra, the method we employ is much faster 
and, in addition, it is easier to implement. 

We construct aggregates with various fractal dimensions using a sequential tunable 
particle-cluster aggregation method developed by Filippov et al. [25]. A fractal 
aggregate composed of homogeneous spheres obeys the so-called scaling law [25] 



Here N is the number of constituents, each with radius a; kf is the fractal prefactor; 
Df is the fractal dimension, and R g is the radius of gyration defined by 



where r, is the position of the ith constituent. The value of the fractal dimension 
can in theory vary between the two extremes Df = 1 (a thin, straight chain of 
particles) and Df = 3 (a homogeneous sphere). The aggregation method we em- 
ploy ensures that with every particle that is added to the aggregate the scaling law 
(Eq. 26) is fulfilled exactly [25]. For all aggregates we choose the fractal prefactor 
kf = 2. Different values of kf result in slightly different fractal aggregate shapes. 
By comparing the results obtained for various aggregates using different values of 
the fractal prefactor we find that the results are not very sensitive to the exact value 
chosen. 




(26) 




(27) 




(28) 
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Fig. 4. Same as Fig. 3 but for fractal aggregates. From left to right the fractal dimension 
Df = 2.8, 2.4, 1.8 and 1.2, respectively. For all four fractal dimensions we considered an 
ensemble of 10 particles with the same shape parameters but different seeds of the random 
number generator. This results in 10 fractal aggregates with the same fractal dimension but 
different locations of the constituents. By averaging over these aggregates we reduce the 
effects of random variations of the particle shapes on the final results. 

The structure of aggregates depends strongly on the precise mechanism by and en- 
vironment in which these aggregates form. The shape of aggregates formed in high 
density environments, for example small gold particles in aqueous suspension, is 
determined by the random-walk nature of the different particles upon approach. 
This process can lead to aggregates with a fractal dimension of about 1.7 [26]. 
In low density environments like protoplanetary disks and molecular clouds in as- 
trophysics, particles approach each other on straight lines before collisions. If such 
clusters grow by addition of single grains to an aggregate ("Ballistic particle cluster 
aggregation", BPCA), this process leads to aggregates with a fractal dimension of 
exactly 3 in the limit of large aggregates [27] , and slightly less than three for smaller 
aggregates. If on the other hand, the growth is approximately mono-disperse, i.e. 
each aggregate grows by collisions with aggregates of similar size ("Ballistic clus- 
ter cluster aggregation", BCCA), the expected fractal dimension is between 1.8 
and 2.1 [28]. Almost linear aggregates (Df « 1) may be formed from particles with 
electric or magnetic dipoles in external fields [29,30] 
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Some typical examples of fractal aggregates are shown above the two panels of 
Fig. 4. The top panel of Fig. 4 shows the ensemble and orientation averaged form- 
factor distributions of ensembles of fractal aggregates with different fractal dimen- 
sions. Each ensemble consists of fractal aggregates with the same fractal dimension 
and fractal prefactor, but with different values of the seed of the random number 
generator. Also shown in Fig. 4 (bottom panel) is the average polarizability per unit 
volume as a function of the real and imaginary part of the refractive index of the 
particles. 

It is apparent from Fig. 4 that the imaginary part of the polarizability per unit vol- 
ume of compact aggregates (Df = 2.8) resembles that of compact homogeneous 
spheres (see the left panel of Fig. 2). When the fractal dimension is decreased, the 
resonance at m = ?' V2~ disappears. Note the strong similarity between the bottom 
panels of Figs. 3 and 4. Although the Gaussian random spheres are very different 
from the fractal aggregates, the absorption properties for both classes of particles 
are quite similar. 



4 Discussion 

4.1 Comparison with distributions of simple shapes 

In order to model irregularly shaped particles, several distributions of simple shapes 
are presented in the literature. We consider here the widely used Continuous Distri- 
bution of Ellipsoids (CDE [13]) and the Distribution of Hollow Spheres (DHS [7]). 
For these distributions we compute the corresponding form-factor distributions and 
compare these to the form-factor distributions of irregularly shaped particles. 

The CDE averages over all possible tri-axial ellipsoids with a specific weighting 
function (see [7]). The sum of the form-factors of a tri-axial ellipsoid averaged 
over all orientations has to be unity. Therefore, there are only two independent 
form-factors, L x and L 2 and one form-factor 1 -L x -L 2 . Since all form-factors have 
to be positive, another constraint is that L 2 < 1 - L x . The CDE averages over all 
possible ellipsoidal shapes by taking the distribution function P(Li,L 2 ) = 2. The 
average polarizability then becomes 




dL 2 (a(L 1 ) + a(L 2 ) + a(l - L x - L 2 )) 



o 



o 




(29) 



o 
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Thus, the form-factor distribution according to the CDE is 



P(L) = 2(1 - L). (30) 

In the DHS we average over the volume fraction, /, occupied by the central vacuum 
inclusion of a hollow sphere, while keeping the material volume of the particle 
constant. This implies that the outer radius of such a particle increases when / 
increases. We consider the shape of a homogeneous hollow sphere to be different 
from that of a homogeneous sphere. As can be verified from Eq. (9) by substitution, 
the form-factor distribution of a single hollow sphere displays two features located 
at 



U= l -- l -j&r^, (3D 
1 1 



L 2 = - + -,JSfTl, (32) 
which have respective strengths 



w x = \ + —^=, (33) 

w 2 = \ p ■ (34) 

From these equations it is already clear that the part 1/3 < L < 2/3 in the form- 
factor distribution is not covered when using hollow spheres. From the above equa- 
tions, the form-factor distribution for the DHS can be obtained and yields 



||2-3L|, 0<L<| and § < L < 1, 

(35) 

0, \<L<\. 



Fig. 5 shows the form-factor distributions of CDE (Eq. 30) and DHS (Eq. 35) to- 
gether with those of an ensemble of Gaussian random spheres with cr = 0.7 in ran- 
dom orientation and an ensemble of fractal aggregates with fractal dimension 1.8 
in random orientation. We can see that the form-factor distribution for the CDE and 
the Gaussian random spheres are in many respects very similar. This might explain 
the success of the CDE as a statistical representation of the absorption properties 
of irregularly shaped particles. We can also see from Fig. 5 that by using hollow 
spheres we do not cover the range of form-factors 1/3 < L < 2/3 (cf. Eq. 35). 
Furthermore, the DHS shows a strong contribution of L > 0.8 which is absent in 
the form-factor distributions of realistically shaped particles. 
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Fig. 5. The form-factor distribution of various particle shape distributions. The solid black 
line represents the distribution for an ensemble of randomly oriented Gaussian Random 
Spheres with cr = 0.7, the solid gray line represents an ensemble of randomly oriented 
fractal aggregates with fractal dimension Df - 1.8, the dotted line represents the Contin- 
uous Distribution of Ellipsoids (CDE [13]) and the dashed line refers to the Distribution 
of Hollow Spheres (DHS [7]). We chose cr = 0.7 for the Gaussian random spheres since 
more modest random spheres (cr = 0.1 or 0.3) still show too much resonance behavior (see 
Fig. 3). The value of the fractal dimension Df - 1.8 is chosen because this is approximately 
expected from cluster-cluster aggregation [28] and it represents a reasonably fluffy aggre- 
gate. Note the remarkable similarity between the Gaussian Random Spheres with cr = 0.7 
and the CDE form-factor distribution. 

4.2 Absorption spectra 



If we consider the wavelength dependence of the refractive index of the particle 
material, we can construct the absorption spectrum of the particle, i.e. its absorption 
cross section as a function of wavelength, using the form-factor distribution. In 
general, this absorption spectrum will display resonances, the exact positions and 
strengths of which are sensitive to the shape of the particles. As an example we 
consider the resonance of a single Lorentz oscillator. Its refractive index is given 
by (see e.g. [13]) 

fto 2 

m 2 = ml + — . (36) 

o>q - or - lyto 



In this equation m > 1 is the real valued refractive index for to —> oo, / is the 
oscillator strength of the feature, to p is the plasma frequency, to is the frequency of 
incident radiation and y is a damping factor. If to is given in wavenumbers, we can 
express Eq. (36) in wavelengths by using to = A' 1 . It can be shown that for a single 
form-factor, L, the absorption spectrum due to a single resonance given by Eq. (36) 
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displays a single feature at position to max given by [13] 



2 2 , 

W max - w + 



(37) 



1 +L(ml- 1) ' 



with a typical width y. From this equation we can see that u> max is an increasing 
function of L. This means that low values of L give a red feature (shifted towards 
long wavelengths, low frequencies), while high values of L result in a blue feature 
(shifted towards short wavelengths, high frequencies). Eq. (37) shows that the posi- 
tion of a resonance, co max , caused by a dust grain in the Rayleigh domain must obey 
the inequality 



The maximum value of the imaginary part of the polarizability for a given Lorentz 
oscillator is also determined by the form-factor L. Furthermore, 



Thus, a low value of L will result in a high value of the imaginary part of the 
polarizability. Combining this with the fact that low values of L give a red feature, 
and that all of the form-factor distributions we computed give more weight to the 
lower values of L than to the higher values, we conclude that, in general, particle 
shape effects tend to a broadening and a red-shift of the spectral features. This is in 
agreement with the findings presented in ref. [7]. 



4.3 The statistical approach 

The method described in this paper provides a strong argument in favor of the 
main assumption of the statistical approach, namely that the optical properties of an 
ensemble of irregularly shaped particles can be represented in a statistical sense by 
the average properties of an ensemble of particles with the same composition, but 
with simple shapes. We have proved that for the absorption properties of particles 
in the Rayleigh domain, the statistical approach has an analytical basis given by the 
form-factor distribution. This implies that the average absorption cross section of an 
ensemble of arbitrarily shaped and arbitrarily oriented particles can be represented 
by a shape distribution of spheroidal particles with the same composition and in a 
fixed orientation. 



a> 2 < a> 2 max < oj + 



(38) 



1 + (m% - 1) " 




(39) 
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We note that the form-factor distribution only provides the absorption cross sec- 
tions, and not the scattering properties of small particles. It is in general not possi- 
ble to find a distribution of spheroidal particles that gives both the absorption and 
the scattering cross section of an ensemble of arbitrarily shaped particles for every 
value of the refractive index. For an ensemble of ellipsoidal particles there is a re- 
lation between the scattering and the absorption cross section which is independent 
of the shape of the ellipsoids and is given by [7] 

Cra . c, b , . (40) 

on Im(m 2 ) 



From Eqs. (21) and (22) it is clear that for arbitrarily shaped particles such a simple 
relation which is independent of the particle shape can only be found in very spe- 
cific cases. This implies that in general there is no shape distribution of spheroids 
that can provide both the absorption and the scattering cross sections of an arbitrar- 
ily shaped particle for every value of the refractive index. 



4.4 Extrapolation to larger particles 



The form factor distribution can be used to construct for a given ensemble of ir- 
regularly shaped particles a shape distribution of simply shaped particles that has 
the same absorption properties in the Rayleigh domain. This distribution of simply 
shaped particles can be used for computations outside the Rayleigh domain. Com- 
bining the form-factor distribution of an irregularly shaped particle, V (L), with the 
equations for the form-factor of a spheroid with the field applied along the rotation 
axis, L(a/b) (see Eq. 8), the shape distribution of spheroids corresponding to this 
irregular particle shape is given by 



P(a/b) = P(L(a/b)) 



dL(a/b) 



d(a/b) 



(41) 



Computations for spheroids in a fixed orientation can be computed relatively fast 
for large particles using for example the Separation of Variables Method (SVM 
[31]) or the T-matrix method [4]. Thus, Eq. (41) provides a well founded choice for 
a distribution of spheroidal particles that can be employed for the computation of 
absorption properties of arbitrarily shaped particles with sizes outside the Rayleigh 
domain. However, it needs to be proven that this approach provides the correct 
optical properties in all cases. 
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5 Conclusions 



We have presented an easy to use method to compute the absorption and scattering 
properties of small particles with arbitrary shape, structure, orientation and com- 
position. The method is based on a solution of the DDA equations in the Rayleigh 
domain. For a given geometrical shape of the particles, the solution has to be com- 
puted only once to obtain the absorption and scattering properties for arbitrary val- 
ues of the refractive index. This provides a significant speedup of the computations 
in cases where calculations for many values of the refractive index have to be done. 
For example, it allows for fast computations of absorption spectra of arbitrarily 
shaped and arbitrarily oriented particles in the Rayleigh domain. Other practical ap- 
plications are, for example, interpretation of atmospheric and radar measurements 
as well as calculations of the scattering matrix of small particles as a function of 
the scattering angle. 

The method we use can be employed for calculating a form-factor distribution. 
This distribution uniquely determines the absorption properties of a particle or an 
ensemble of particles. Using the form-factor distribution, we have studied the basis 
of the statistical approach in the Rayleigh domain. Various shape distributions of 
particles can then exhibit the same absorption properties. Two ensembles of dust 
grains with the same composition and the same average form-factor distribution, 
but with different shape distributions have exactly the same absorption properties. 
Therefore, it is not possible to distinguish between these two ensembles from their 
absorption properties alone. Furthermore, for a given form-factor distribution it is 
trivial to obtain a distribution of spheroidal shapes with exactly the same average 
form-factor distribution. This provides a strong argument in favor of the fundamen- 
tal assumption of the statistical approach, i.e. that the average absorption and scat- 
tering properties of an ensemble of irregularly shaped particles can be represented 
by an ensemble of simple shapes. 

The form-factor distribution uniquely determines the absorption properties of ar- 
bitrarily shaped particles allowing for the construction of a shape distribution of 
spheroids with the same composition and a fixed orientation with the same av- 
erage absorption cross section. It is, however, in general not possible to obtain a 
shape distribution of spheroidal particles with both the same average absorption 
and the same average scattering cross section. It needs to be investigated what the 
error on the computed scattering properties is when using the shape distribution of 
spheroids. 

The conclusions presented in this paper are valid for particles in the Rayleigh do- 
main. The form-factor distribution of an ensemble of irregularly shaped particles 
can be used to construct a shape distribution of spheroidal particles in a fixed ori- 
entation with exactly the same absorption properties in the Rayleigh domain. Since 
the absorption properties of spheroidal particles in a fixed orientation can be ob- 
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tained relatively easy also for larger particles (i.e. outside the Rayleigh domain), 
this shape distribution provides a well founded choice for studying combined par- 
ticle size and shape effects. 
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